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Nuclear spin-lattice relaxation rate (1/Ti) measurements in optimally doped Ba(Fei_i,Rlia;)2As2 
superconductor evidence the onset of low-frequency dynamics associated with the vortex motions. 
The slowing down of the dynamics yields a peak in l/Ti, which has been observed to progressively 
shift to lower temperatures as the magnetic field is increased from 3 T to 15 T. We interpreted 
these results in the light of a phase transition between a liquid vortex state and a glassy vortex 
ground-state. The Arrhenius law, employed to describe the vortex dynamics in the cuprates, fails 
in the case of the iron-pnitictides, whereas a Vogel-Tammann-Fulcher law gives more meaningful 
results, together with the dynamical scaling approach. The glassy phase transition occurs below Tc, 
very close to the Kosterlitz-Thouless transition temperature of the vortex lattice. We further show 
that the anisotropy in the vortex dynamics causes also a dependence of the spin-spin relaxation rate 
on the magnetic field orientation. 



PACS numbers: 74.25.nj, 74.25. Wx,74.25.Uv 

I. INTRODUCTION 

The discovery of iron-based superconductors ^ has re- 
newed the interest in the vortex physics. Thanks to their 
layered structure, the pnictides are good candidates for 
studying the dynamical properties of the flux line lattice. 
In particular the 122 family of these new superconduc- 
tors can be grown in large high-quality single crystals?^ 
suitable for nuclear magnetic resonance (NMR) experi- 
ments. On the other hand, since iron-pnictides are not 
as anisotropic as the cuprates or as the organic super- 
conductors, one may wonder if the liquid vortex phase 
is extended enough to be investigated, and whether a 
transition between a liquid and a solid phase or rather a 
glassy-phase would occur in these novel materials. 

More than 40 years ago it was argued that the material 
imperfections destroy the crystalline long-range order of 
the Abrikosov vortex lattice."^ The quenched disorder is 
responsible for the pinning of the vortices and for the 
contextual disruption of the lattice structure.'* After this 
theoretical prediction it was experimentally found^^— that 
in the cuprates the flux-lattice melts in a glassy state oc- 
curring well below the transition temperature Tc- More 
recent ac-susceptibility measurements seem to evidence 
a glassy state in the 122 family of iron-pnictides, too^. 
Some open questions concern the possibility of observing 
a glassy phase transition in a 3D array or vortices, as well 
as the nature of such a transition: whether it is a true 
thermodynamic transition with a defined order parame- 
ter, or a cross-over between two regimes Another issue 
concerns the effects of the weak anisotropy on the geom- 
etry of the vortices: whether they have to be thought as 
2D vortices (pancakes) or rather proper flux tubes, and 
how the phase transition between 2D and 3D lattice can 



be driven. 

The investigation of the vortex motion is important 
both for the preservation of the dissipation-less proper- 
ties of superconductor^ and for fundamental physical as- 
pects involved. The phenomena hidden in vortex physics 
are extremely rich, and can respond to many experimen- 
tal techniques. Apart from the pinning mechanism, pro- 
vided by lattice imperfections or vacancies, there is an- 
other source of pinning, the intrinsic pinning, due to the 
anisotropy of these systems. In case of a layered struc- 
ture, when the applied field is in the c direction, the vor- 
tices are not affected by the anisotropy in the ab plane, 
but when the field is applied in the ab plane, the vortex 
cores are located between the layers, in order to mini- 
mize their energy i^i In the absence of pinning centers, 
vortex motion within the layers is allowed, but motion 
perpendicular to the layers is hindered. 

Recent nuclear magnetic resonance (NMR)i^ii^ and 
muon spin rotation (/iSR)^ measurements on iron pnic- 
tides have studied the dynamics of the flux lines lattice, 
but to the author's knowledge, no systematic study of 
both the temperature and of the magnetic field effect 
has been done, so far. 

This paper extends our previous work (Ref. IT^ where 
the NMR spin-lattice relaxation rate was measured at 
7 T and 3 T, in a Rh-doped member of the 122 family 
of iron-based superconductors, near the optimally doped 
regime (Ba(Fei_j;Rh3;)2As2, with X'-^ 0.07). We show 
how the spin-lattice relaxation rate of the ""^As nuclei 
changes with temperature and field, and propose an in- 
terpretation based on a 3D flux line lattice (FLL) model, 
where the thermal excitations induce fluctuations per- 
pendicular to the external magnetic field. We also notice 
that the Arrenhius law fails in describing the correlation 
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time of the vortex motion and that the Vogel-Tammann- 
Fulcher (VTF) law should be rather employed. This sug- 
gests the occurrence of a glassy phase transition. We 
discuss further on the I/T2 relaxation rate measured at 
high field, and its dependence on the field orientation. 



II. TECHNICAL ASPECTS AND 
EXPERIMENTAL RESULTS 

The ^^As NMR measurements were performed on a flat 
0.8 X 5 X 7 mm'^ single crystal of Ba(Feo.93Rho.o7)2As2 
(BaPeRAs), with the c axis along the shortest side. The 
details about the sample are presented in Ref. [T3. The 
earlier experiments performed on the same sample ev- 
idenced a peak in the spin-lattice relaxation rate, be- 
low the transition temperature Tc, whose position and 
amplitude strongly depend on the magnetic field {H) 
amplitude. The lower field experiments {H = 3 T, 7 
T) are here completed with higher field ones. Exper- 
iments at = 8.5 T were performed in Pavia, while 
the 11 T and 15 T experiments were carried out at the 
LNCMI-Grenoble. The results show remarkable similar- 
ities with an earlier study on oriented powders of the 
cuprate YBa2Cu408 (YBC0124), with a critical tem- 
perature ~ 81 K, and a higher anisotropy ratio 7.^^ 

The NMR measurements were performed by using 
standard radio frequency pulse sequences. The spin- 
lattice relaxation time Ti was measured by means of 
a saturation recovery pulse sequence and the recovery 
of the nuclear magnetization m{t) was well fitted by a 
slightly stretched double exponential law, corresponding 
to a magnetic relaxation mechanism of a nuclear spin 
/ = 3/2, for the central (—1/2, 1/2) transition: 

l-m{t)/ma = 0.1exp(-V?'i)'' + 0-9exp(-6i/Ti)'', (1) 

with r -> 1. Pig. [T] shows the relaxation rate T^^ at dif- 
ferent fields, with H || c, where a peak is observed below 
Tc- We immediately notice that, by increasing the field, 
the peak in l/Ti is rapidly reduced above 7 T, so that 
only a small kink remains at 11 T and 15 T. As shown in 
Ref. [l^ the peaks are visible just in H || c geometry (Fig. 
1). The contribution of the vortex dynamics to the spin- 
lattice relaxation rate \/Tivl is found by subtracting 
from the raw data the electronic contribution, estimated 
by extrapolating the data out of the peak region. This 
procedure is in practice equivalent to subtracting from 
the raw data the ones taken for H _L c, properly rescaled 
by the hyperfine factors, since for that orientation the 
nuclear spin-lattice relaxation is only given by the elec- 
tronic spin fluctuations (Pig. [2]). Namely, one has 
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where the first term is the vortex-lattice contribution, 
while the second one is due to electron spin fluctuations. 
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FIG. 1. The spin-lattice relaxation rate at different magnetic 
field values: 3 T (blue squares), 7 T (purple triangles), 8.5 T 
(green circles), 11 T (yellow stars), and 15 T (pink pentagon). 
All data refer to H || c orientation. The arrows show the 
temperature of the detuning of the NMR probe resonating 
circuit. 
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FIG. 2. The spin-lattice relaxation rate measured at 7 T, for 
H II c (black diamonds) and for H _L c (blue circles), rescaled 
by the hyperfine coupling ratio (see Refill) . The red dotted 
line is the interpolation of the data out of the peak region. 
The peak is visible just in the longitudinal field geometry, 
while in the other one, just the electronic contribution 
to the nuclear spin relaxation is present. This result, already 
shown in Ref. [12 . has been report here for the sake of clarity. 



Further information on the vortex dynamics can be ob- 
tained from the spin-echo relaxation time (T2) measure- 
ments below Tc, at 15 T. To measure T2 we employed 
a Hahn-echo sequence, by varying the delay time r be- 
tween the two pulses. By dcflning T2 as the time at 
which the echo drops at 1/e of its initial value, we observe 
a pronounced dependence of the relaxation rate on the 
direction of the applied field, with respect to the FeAs 
planes: when H || c the spin-echo decay rate increases 
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with decreasing temperature. On the other hand, when 
H _L c, the rate saturates to a value in agreement with 
the calculation of the dipolar contribution to the relax- 
ation following Van-Vleck summation&i^ (See Fig. [3]). 
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FIG. 3. The figure shows the spin-echo decay rate at 15 T, 
for H ± c (red circles) and H || c (black squares), determined 
by time where the echo decays to 1/e of its initial value. 



Nevertheless, we didn't observe any sudden transition in 
the NMR relaxation times. In addition, below the tem- 
perature of the peak in a marked increase in the 
acoustic ringing was observed at all field values, mostly 
when H || c. Now we will only briefly comment on this, 
since this effect will be the subject of a future work. We 
notice that a strong ringing was observed in Fe-doped 
2H-NbSe2 as well,^^ in the proximity of the spin-glass 
phase transition. We argue that if a glass-phase tran- 
sition occurred the vortices, that in the liquid phase are 
free to move reversibly, would transfer their elastic energy 
to the lattice (via the pinning mechanism) . The compar- 
ison with the spin-glass system can be understood if we 
assume that a coupling between the spin fluctuations and 
the phonons is also playing a role. 

Furthermore, when studying layered superconductors, 
it is necessary to deal with the dimensionality crossover 
that might occur in the FLL. The crossover field Bcr 
was calculated by L. I. Glazman and A. E. Koshclcv. 
In the limit Ac < ^ab/^^ with Ac and Xat the London 
penetration depth for H J- c and H \\ c, respectively, 
namely when the Josephson coupling dominates over the 
magnetic coupling, one can writ o^^i^° 
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III. DISCUSSION 

Let us start our analysis with some considerations on 
the phase diagram of our compound. Generally, the 
magnetic phase diagrams for type-II superconductors are 
rather complex: in addition to the complete Meissner 
phase, a region with pinned or otherwise non-moving vor- 
tices is predicted at very low temperature, whereas at 
higher temperatures the FLL melts into a liquid phase, 
where vortices are free to move reversibly. The bound- 
ary line between the two phases is the irreversibility line, 
sometimes also called the melting line. Actually the irre- 
versibility and melting do not always coincide, and Ma- 
jer et al.^ showed that in Bi2Sr2CaCu208 (BSCCO) the 
two processes stem from different and unrelated physical 
mechanisms. The irreversibility line defines the boundary 
between irreversible and reversible motions of the vortices 
and should depend on the time scale of the experimen- 
tal technique, as well as on the geometry of the sam- 
ple, while the melting line is independent of the experi- 
ment. Above this line, the vortex cores fluctuations can 
be strong enough to cause the vortex solid or glass to melt 
into a liquid phase. When the melting occurs, according 
to the Lindemann criterion, the root-mean square ampli- 
tude of the vortex displacement < >i/2 exceeds leCL, 
where Ig is the vortex lattice step and cj^ — 0.1 — 0.2 is the 
Lindemann parameterjii In case of a melting transition 
an abrupt change in the NMR quantities, a sharp peak 
or a step in the spin-lattice relaxation rate is expected, 
owing to the first-order nature of the phase transition. 



where s is the inter-layer distance, $o is the quantum 
of magnetic flux, ^ is the coherence length. Since here 
7 ~ 2 we calculated a critical fleld of about 90 T, which 
implies that the FLL is in the 3D correlated regime, for 
all the applied fields. At such low fields the usual picture 
of "flux lines" can be restored. 

Another important effect is due to the temperature. 
In 1991 R. J. Clem^^ predicted a divergence of the root 
mean square amplitude of the fluctuations, for an isolated 
vortex line, given by a stack of 2D pancakes, in absence 
of the inter-layer Josephson coupling, at the temperature 
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When the vortices density is large enough that an array 
of flux lines enters into the system, the thermal energy 
can still decouple the vortices in different layers. The new 
crossover temperature T^^ depends on the field, but in the 
H ^ limit, the equivalence T^^ = Tcr is restored. Tcr 
was identified by Clem^ii^ with the Kosterlitz-Thouless 
transition temperature. This crossover was studied in 
the cuprates by employing ^SR,^'^ and NMR^i. We cal- 
culated that for our compound Tcr ^ 17 K, by taking into 
account the temperature dependence of London penetra- 
tion depth A. 

As far as the peak in the 1 /Ti is concerned, we started 
our analysis with a BPP-like relation for the spectral den- 
sity of the spin fluctuations.^^ As shown in Ref. 1/Ti 
due to vortex motions can be written in the forroi^: 
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where ^"^7 is the As gyromagnetic ratio and < >3d 
the root mean-squared ampUtude of the local field fluc- 
tuations when the vortex lattice is correlated along the 
c-axisJ^ le is the inter-vortex spacing, Tc is the correla- 
tion time, which depends on the dynamics involved in the 
fluctuations, and lol is the Larmor frequency. We firstly 
assumed that Tc had an Arrenhius-like form, as claimed 
in YBa2Cu40g^'^ and HgBa2Cu04_5^^ cuprates, where 
a collective activated motion is thought to describe the 
dynamics acting in the FLL: 
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By using the above equation to analyze BaFeRAs data, 
we obtained non-physical results for the intra-valley cor- 
relation time To (10~^^ s). An anomalous value of tq was 
firstly observed in our previous work,-'^^ 

We then considered the VTF law (Fig. H]), employed 
to describe the divergence of the relaxation time in a 
glass-forming system: 



Tc{T,H) = Toexp 
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where U^ff refers an effective energy barrier and Tq is 
the temperature where the divergence occursi^"— This 
interpretation was suggested by the theoretical prediction 
of Fisher et al.'^ In this case a divergence of the viscosity 
7] can be related to the divergence of the correlation time 
77 = rGoo, where Goo is the high frequency shear mod- 
ulus. With this correction the fit gives more meaningful 
values for the intra-valley correlation times tq . 

To test the validity of this approach we employed the 
VTF law to fit the peak in 1/Ti in the cuprates (Fig. [5]) 
and we obtained quite nice results for the lowest fields. 
However, in this case also the Arrhenius law gave good re- 
sults, so it is not easy to distinguish between the two sce- 
narios. We recall that the Arrhenius law is employed to 
describe the divergence of the viscosity in glass-forming 
systems, like silica, which are called "strong" systems, 
whereas the VTF law is rather employed to describe 
the "fragile" systems;^ The terminology "fragile" and 
"strong" refers to the extent in which the temperature 
dependence of the viscosity deviates from Arrhenius be- 
havior: the fragility m of glass-forming systems refiects 
the slope (i.e. logarithmic derivative) of the viscosity of 
a material with temperature, as it approaches the glass 
transition temperature: 
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where r is the relaxation time and the derivative is eval- 
uated at T = Tg. We further comment on the width of 
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the peaks in the cuprates by observing that the sample 
was an oriented powder, instead of single crystal, where 
a distribution of correlation times might have broadened 
the peak. 
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FIG. 4. The figure shows the spin-lattice relaxation rates, at 3 
T (green circles), 7 T (blue squares), 8.5 T (black diamonds), 
for H II c. The electronic contribution has been subtracted 
from the raw data. The red lines are the best fits, according 
to Eq.s[5]and[7l 

A second approach to describe the glass-phase tran- 
sition is the dynamical scaling. The vortex-glass corre- 
lation length ^DSi which is the length scale over which 
the phase of the pair wave function is correlated, and the 
characteristic relaxation time r^s both diverge at the 
glass transition: ^'''^^ 
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where Tg is the glass-phase transition temperature and 
P the corresponding critical exponent, and e the reduced 
temperature. In Table I we present the fitting values, 
with their error bars, for both the VTF law and the 
power-law fitting parameters. We notice that the glass 
phase transition temperature Tg is in good agreement 
with the instability temperature Tqi extracted by the 
VTF law. In addition, the /3 exponent agrees with the 
predicted universality class for 3D systems, where it is 
expected to be 6 ± \P 

If we collect these results in a phase diagram (Fig. [6]) , 
we notice that the "instability-line" intercepts the tem- 
peratures axis near T^r ^ 17 K, which marks a crossover 
between a regime where the vortices are correlated and a 
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FIG. 5. The figure shows the spin-lattice relaxation rates 
data, at 5.9 T (green diamonds), 9.4 T (blue squares) and 14.9 
T (black circles) in YBC0124, when H || c. The electronic 
contribution has been subtracted from the raw data. The 
data in the figure are taken fi-om Ref. [3^ , whereas the red 
lines correspond to a fit using a correlation time described by 
a VTF law. For the case of 5.9 T and 9.4 T the fit quality 
is satisfactory, while this is not the case for the highest field 
(14.9 T) measurement (see text). 

high temperature one, where they are uncorrelated. Fur- 
thermore the fitting results and the peak widths appear 
to be in agreement with a liquid-glassy phase transition 
rather than with a first order transition. We stress that 
the line connecting the position of the peak in the 1/Ti 
doesn't mark any phase boundary, but on the other hand 
it marks the temperature at which the correlation time 
matches the inverse of the Larmor frequency, leading to 
an enhancement of the spectral density of the transverse 
field fiuctuations. 



H (T) 


Ueff(K) 


To(K) 


Tg(K) 


P 


3 T 


22.4 ± 12.4 


16 ± 1.6 


16.5 ±0.3 


5.9 ± 0.4 


7 T 


19.7 ± 3 


15±0.4 


15.5± 0.5 


6.3±0.9 


8.5 T 


16.9±4.9 


14.5 ±0.6 


15±0.2 


5.8±03 



TABLE I. The table presents the fitting values, for the VFT 
law and the dynamical scaling approach. In both cases tq 
is about 10~^^ s, as expected for the correlation time for the 
vortex intra- valley motion.— Similar time values have been re- 
cently obtained by ac-susceptibility measurement on the Co- 
doped member of the same family of iron-pnictides^. 



Let us now comment on the height of the peaks: both 
the cuprate and the pnictide show an increasing ampli- 
tude of the l/Ti peak around 35 K, for YBC0124, and 
around 17 K, for BaFeRAs. This particular tempera- 
ture coincides with Tcr for both cases. This may ex- 
plain the behavior of such peaks with the magnetic field: 
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FIG. 6. The equilibrium phase diagram (j = 0) for the vortex 
lattice. The black squares show the upper critical field Hc2 
line, measured by SQUID magnetometry. The red solid line 
is a guide to the eye. The black circles refer to the peak in the 
spin-lattice relaxation rate. The green dotted line is an guide 
to the eye which does not mark any phase boundary (see text). 
The blue half full circles mark the instability temperature To, 
in agreement with the VTF law and the blue arrow marks Tcr- 
The blue solid line is a guide to the eye. The dimensional 
crossover field is marked with an arrow, because it is 

expected at very high fields. 



by applying higher fields, the peaks move towards lower 
temperatures, so that the correlation time, and accord- 
ingly the spin-lattice relaxation rate, are strongly affected 
by the VTF law, via equation [S] We further observe 
that the glassy-phase transition temperature at zero field, 
To(0) ^ Tg[Q), is expected below Tc, as observed in the 
cuprates.*'^ Remarkably, this glass-transition tempera- 
ture is found very close to the Kosterlitz-Thouless tem- 
perature. The same result was also found in the cuprates, 
where ~ 35 K. 



IV. CONCLUSIONS 

This work completes a previous NMR study on a su- 
perconductor belonging to the 122 family of iron-based 
superconductors. We investigated the fluctuations of lo- 
cal fields and we interpreted them in terms of the vortex 
dynamics. We have studied how the magnetic field in- 
tensity affects the peak in the spin-lattice relaxation rate 
1/Ti, and we found that this behavior is consistent with 
a liquid-glassy phase transition occurring in the vortex 
lattice. Our results are compared with the ones for the 
YBC0124 cuprate, where similar peaks were found. In- 
formation on the scaling exponents for the dynamics and 
on the correlation time was obtained and shown to be 
consistent with a 3D vortex lattice model. 
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